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Abstract. Pulsar glitches are a potential source of gravitational waves for 
O , current and future interferometric gravitational wave detectors. Some pulsar 

C^i glitch events were observed by radio and X-ray telescopes during the fifth LIGO 

science run. It is expected that glitches from these same pulsars should also 
be seen in the future. We carried out Monte Carlo simulations to estimate the 
sensitivity of possible gravitational wave signals associated with a pulsar glitch 
using a coherent network analysis method. We show the detection efficiency and 
evaluate the reconstruction accuracy of gravitational waveforms using a matched 
^ ■ filter analysis on the estimated gravitational waveforms from the coherent analysis 

*yv ' algorithm. 
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PACS numbers: 04.80.Nn, 07.05.Kf, 95.55.Ym, 95.85.Nv, 95.85.Sz 



1. Introduction 



Pulsar glitches are a sudden spin-up of a neutron star followed by a long term 
relaxation period. The increase in rotational energy during such a pulsar glitch is 
^^ . typically 10^^ ergs. The time scale of the sudden spin-up is observationally limited 

5h ' to within two minutes [T. The mechanism is understood in terms of spin-down 

induced quakes in the neutron star's crust or the transfer of angular momentum from a 
superfluid in the inner crust. Some models have been proposed to explain the physical 
mechanism of the glitches [21 El HI O [6] , however, the relevance of these models is 
still not certain. Though one needs more understanding of neutron star superfluidity, 
the drastic energy transfer should cause excitation of acoustic modes, resulting in the 
emission of gravitational waves. From the observations of pulsar glitches, one could 
infer not only the mass and radius of the stars, but also the superfluid nature, which 
could give us insight into the mechanism of the pulsar glitches. According to [7], the 
detection of gravitational waves caused by the superfluid modes excited by glitches 
gives us a good chance to understand the structure of the neutron star. 

From theoretical and observational studies, the time intervals from one glitch to 
the next from particular pulsars can be predicted [H] ^. The paper by Middledich 
et al. reports on more than seven years of monitoring PSR 0537-6910, the 16 ms 
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pulsar in the Large Magellanic Cloud by using data acquired with the Rossi X-ray 
Timing Explorer. During the monitoring of this pulsar, 23 glitches were observed. 
From these observations, they found that the time interval from one glitch shows a 
strong linear correlation with the amplitude of the first glitch, with a slope of about 
6.5 days per /zHz, and this interval can be predicted to within a few days. This 
predictability provides a novel opportunity for gravitational wave detection during the 
future LIGO J^ and VIRGO [H] science runs. The prediction of the next glitch from 
an observation of a recent glitch event enables us to adjust the observation schedule for 
world-wide gravitational wave detectors to catch the event. This predictability differs 
from the observations of other gravitational wave sources which allow only accidental 
observations. Two glitches from the pulsar J0537-6910 (Aug. 4, 2006 and Nov. 25, 
2005) and one glitch from the pulsar B0833-45 (Aug. 12, 2006) have been observed by 
radio and X-ray telescopes during the 5th LIGO science run (S5) [9l|T2]. According 
to figure 2 in [5] , it is quite possible that several glitch events will be observed during 
the next LIGO science run. 

We propose to look for gravitational wave signals associated with pulsar glitches 
using a coherent network analysis algorithm. Other proposed searches for gravitational 
waves from pulsar glitches are described in |13j . Unlike our coherent method, which 
does not assume any particular gravitational waveform, the method in [13j looks for 
gravitational waves which have f-mode ringdown waveforms by using Bayesian model 
selection. In a coherent network analysis algorithm, one combines the data streams 
from multiple gravitational wave detectors coherently, taking into account the antenna 
patterns, geographical locations of the detectors, and the sky direction to a source. 
This method takes full advantage of the global network of interferometric gravitational 
wave detectors currently in operation, and improves directional searches, resulting 
in enhanced detection efficiency. Mathematically, the extraction of a gravitational 
wave signal is an inverse problem. Due to the degrees of freedom of a matrix of 
antenna patterns, the matrix becomes rank deficient at some detector networks and 
sky locations, resulting in the amplification of variance of the gravitational wave 
parameters. Another way to understand this is that because of the rank deficiency, 
one of the reconstructed polarization waveforms becomes noisy as the noise level is 
dependent on the degree of the rank deficiency. To solve this ill-posed problem in 
the coherent network analysis, we use Tikhonov regularization [HI [15]. Besides the 
Tikhonov regularization, we added new time domain data conditioning [16j [17j to the 
analysis pipeline, creating a complete stand-alone coherent network analysis pipeline 
called RIDGE [TB]. In this paper, we demonstrate the search for gravitational waves 
associated with pulsar glitches using RIDGE and show the detection efficiency using 
simulated detector noise. 

2. Overview of the RIDGE pipeline 

RIDGE is a coherent network analysis method, which is described in detail in |18j . 
The pipeline consists of two main components, data conditioning and generation of 
detection statistics. 

The aim of the data conditioning is to whiten the data to remove frequency 
dependence and any instrumental artifacts. In RIDGE, the data is whitened by 
estimating the noise floor using a running median [17j . By using a running median, 
the estimation of the noise floor can avoid being affected by any large outlier caused by 
a strong sinusoidal signal. The whitening filters are implemented using digital finite 
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impulse response filters with the transfer function |T(/)| — l/-y/S'(/), where S{f) 
is the power spectral density (PSD) of detector noise as a function of frequency /. 
The filter coefficients are obtained by using the PSD obtained from a user-specified 
training data segment. The digital filter is then applied to a longer stretch of data. 

Narrow-band noise artifacts known as lines are a typical feature of the output of 
interferometric gravitational wave detectors. Lines originate from effects associated 
with the functioning of the detectors such as mirror/suspension resonant modes, power 
line interference during operation, calibration lines, and effects from environmental 
sources such as vacuum pumps [TO]. These lines dominate over a specific band, which 
reduces the signal-to-noise ratio and the accuracy in recovering gravitational wave 
signals. In RIDGE, the whitening step above is followed by a line estimation and 
removal step using a method called median based line tracker (MBLT, for short) and is 
described in [16]. Essentially, the method consists of estimating the amplitude and 
phase modulation of a line feature at a given carrier frequency. The remarkable feature 
of this method is that the line removal does not affect the signals in any significant way. 
This is an inbuilt feature of MBLT which uses the running median for the estimation 
of the line amplitude and phase functions. Transient signals appear as outliers in the 
amplitude/phase time series and are rejected by the running median estimate. 

The resulting conditioned data is passed on to the next step which consists of 
generation of a detection statistic. The basic algorithm implemented for the coherent 
network analysis in RIDGE is Tikhonov regularized maximum likelihood, described in 
[15] . Recent studies [2Ql|T5l|2T] show that the inverse problem of the response matrix of 
a detector network becomes an ill-posed one and the resulting variance of the solution 
becomes large. This comes from the rank deficiency of the detector response matrix. 
The strength of the rank deficiency depends on the sky location, and therefore, plus 
or cross-polarized gravitational wave signals from some directions on the sky become 
too noisy. In RIDGE, we reduce this ill-posed problem by using the Tikhonov regulator, 
which is a function of the sky location. The input to the algorithm is a set of equal 
length, conditioned data segments from the detectors in a given network. The output, 
for a given sky location 9 and </>, is the value of the likelihood of the data maximized 
over all possible /i+ and /ix waveforms with durations less than or equal to the data 
segments. The maximum likelihood values are obtained as a function of 9 and (j) - 
this two dimensional output called a skymap. Using the entire skymap, we construct 
radial distance which scales the maximum likelihood by the mean location of the same 
quantities in the absence of a signal. Detailed description about the detection statistic 
is given in [TS] . 

3. Simulations and results 

We performed Monte Carlo simulations to estimate the detection efficiency of possible 
gravitational waves from pulsar glitches. The network consisted of the 4 km and 
2 km LIGO Hanford (HI and H2) and LIGO Livingston (LI) interferometers. For 
the detector noise amplitude spectral densities, we used the design sensitivity curves 
for the LIGO detectors as given in [l^ and kept the locations and orientations 
the same as the real detectors. For H2, the sensitivity is \/2 times less than HI. 
Gaussian, stationary noise was generated (~ 2000 sec) by first using 3 independent 
realizations of white noise and then passing them through finite impulse response 
(FIR) filters having transfer functions that approximately match the design curves. 
To simulate instrumental artifacts, we added sinusoids with large amplitudes at 



Searches for gravitational waves from pulsar glitches 




frequency [Hz] 



Figure 1. The normalized sensitivity and root-sum-square of h of the injected 
ringdown signals used in the Monte Carlo simulations. The open circles are the 
injection signals of /irss = 2.1 X lO^^'^Hz" ' with center frequencies indicated in 



(54, 60, 120, 180, 344, 349, 407, 1051) Hz. 

Signals of constant amplitude were added to the simulated noise at regular 
intervals. The injected signals corresponded to a single source located at the 
right ascension (RA) of 8.6 hours and the declination (DEC) of —45.2 degrees, 
which was the location of Vela pulsar. We assumed that the h+{Q, fc,t) = 
Aex-p{~TTfct/Q)sm{2Trfct) and hx{Q, fc,t) — ^exp(— vr/ct/Q) cos(27r/ci), where t is 
time and t > 0, /c is the central frequency, Q the quality factor and ^ is a constant 
value. The signal strength is specified in terms of root-sum-square defined as /irss = 

-1 1/2 

J^^dt {h\{t) + hl{t)) . In this simulation we took /i^ss = 2.1 x lO^^ijjz-i/s^ 

Q = 5.6 and /^ = 54.7, 62.6, 73, 87.6, 109.5, 146, 219, 437.9 Hz. Unfortunately, 
there is no reliable estimation of the amplitude yet. In this paper we use amplitudes to 
check the performance of this method for the ringdown with wide range of fc- Figure [1] 
shows the normalized sensitivity of the simulated detector noise for HI, H2, and LI 
along with the h^ss of the simulated signals. The normalized sensitivity is defined as 
the amplitude spectral density divided by (F^ -I- F^y^^ where F^, F^ are antenna 
patterns for /i^, /ly . Though this paper uses specific ringdown waveforms with specific 
fci Q- value, this method can recover any type of burst signals since this method does 
not assume any particular model for the signal. Figure [5] shows the averaged antenna 
pattern i^av defined as Ff_ + F^. When the Vela pulsar glitch occurred at 21:59 on 
12th August, 2006, the values of the detector antenna function were 0.16 and 0.12 for 
the LIGO Hanford and Livingston detectors respectively. 

The simulated data was generated at the sampling frequency of 16384 Hz and 
then passed through the data conditioning pipeline. Besides downsampling the data 
to 2048 Hz by applying the same anti-aliasing filter to all data streams, the data 
conditioning pipeline applies time domain whitening filters that were trained, as 
described earlier, on the first two seconds of data for each detector (without any 
injected signals). In real data analysis, we use several seconds of data, which is 
stationary, around the event excluding on-source data. On-source data refers to a 
data set in which the signal is likely to be present. In the case of the VELA pulsar 
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Figure 2. The averaged antenna pattern, -Fav, at the direction to the Vela pulsar 
around the time when the glitch occurred at 21:59 on 12th August, 2006 which 
corresponds to 12 hours in x-axis in the figure. The x-axis start with preceding 12 
hours from the glitch event. During the Vela pulsar glitch, the detector response 
is 0.157 and 0.127 for the Hanford and Livingston detectors respectively. 



glitch, which occurred on Aug. 6, 2006, the reported gUtch epoch is ~ 100 seconds. 
The delay due to dispersion of the radio wave signal through the interstellar medium 
is of the order of 100 ms and can be ignored. Therefore, on-source for the pulsar 
glitch is 100 ms of data around the event. Since we don't know the signal model 
for pulsar glitches, it is difficult to set an optimal length. So we use 0.5 seconds for 
one demonstration. Note that in real data analysis, we should use several integration 
lengths or clustering approach to improve the detection efficiency. 

3.1. ROC curve 

Figure [3] shows Receiver Operating Characteristic (ROC) curves for the simulations 
described above. We calculated ROC curves for the 8 injection ringdown signals which 
are characterized by /ij-ss, the central frequency and the Q-value shown above. The 
detection efficiencies at false probability 5 x 10^'^ are 0.99, 1, 1, 1, 1, 0.79, 0.14 from 
high central frequencies. In regard to the computation cost, all the processing of the 
simulation for one of injection signals presented here took 2 hours using a PC with 2 
GHz CPU. The on-source data can be estimated ~ 120 sec, and even if one hour data 
around the events is taken as off-source data, it's still not computationally expensive. 



3.2. Reconstructed waveform 

Figure |4] shows the reconstructed waveforms hj^, hx for the injection ringdown signals 
with fc — 109.5 Hz. The reconstructed waveforms are subject to noise which 
comes from the detectors. It is necessary to de-noise the waveform when we obtain 
astrophysical information from unmodeled gravitational waves, and some methods for 
this purpose are discussed in [53j. 

Here, assuming the signal waveforms are known, we carry out a matched filter 
method to evaluate the accuracy of the reconstruction of ft,+ , hx- Figure [5] shows a 
flow chart of the pipeline for the matched filter used here. The reconstructed /i+, 
hx are whitened by the whitening filter implemented in RIDGE. The filter parameters 
are estimated using the preceding 2 seconds of data, which do not contain signals. 
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Figure 3. Receiver Operating Characteristic. The sensitivity curve is 
the same as Figure [2] hrss of the injection ringdown signals are 2.1 X 
10~^^Hz~^'^ and characterized by the center frequency and Q-value (fc,Q) = 
(62.6, 5.6), (73, 5.6), (87.6, 5.6), (109.5, 5.6), (146, 5.6), (219, 5.6), (437.9, 5.6). 
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Figure 4. Plots showing the reconstructed /i+ (left) and hx (right). The 
corresponding injection ringdown signal is the one with a central frequency of 
109.5 Hz, Q-value of 5.6, and hrss is 2.1 X lO-^lRz-^/^. 
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Figure 5. Block diagram of the matched filter method which is appUed to the 
reconstructed /i_|_, h^. The /i_|_, h^ reconstructed by the RIDGE pipeline are 
passed on to the next step of whitening the reconstructed data streams. The 
filter coefficients are estimated using the preceding data segment, which does not 
contain injection signals. The template signal is also whitened by the filters with 
the same filter coefficients. The matched filter method is applied to the resulting 
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Figure 6. Upper plots show the reconstructed h+ (left), h^ (right) whitened by 
the data conditioning routine in RIDGE. Lower plots show the output of matched 
filtering for reconstructed h^, hx The output is the signal-to- noise ratio (SNR). 



The template s is also whitened by the same whitening filter and filter parameters. 
The whitened data is passed on to the matched filter method which maximizes 

T = x^ s'' /aV s^ s'f, where x^ and s^ are column vectors of the whitened data hi, ^, 
and whitened template respectively, a is the standard deviation of the whitened data. 
T corresponds to the signal-to-noise ratio (SNR). 

Figure [6] is the result. Upper plots show the whitened, reconstructed h+, hx- 
Lower plots show that the SNR of the reconstructed h+, hx are 6.9, 6.3. Though any 
coherent network analysis has bias and we don't correct the bias, this result is not bad 
compared with the ideal SNRs of the injection signal: 10.1 for HI, 6.8 for H2, 8.4 for 
LI. 



4. Conclusion 

Pulsar glitches are a potential source of gravitational waves, and some pulsar glitches 
were observed during S5 by radio and X-ray telescopes. Our work is in progress 
to search for possible gravitational wave signals associated with all the observed 
pulsar glitches. The next glitch from a given pulsar can be predicted from previous 
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observations, which gives us a chance to adjust the observation schedule, so that the 
detectors are on when the ghtch occurs. This predictabihty of the pulsar glitches 
creates a new type of gravitational wave search as compared to other triggered 
gravitational wave searches which only allow serendipitous detections. In this paper, 
we perform Monte Carlo simulations to study the detection efficiency of possible 
gravitational waves triggered by the Vela pulsar glitch on 12th August 2006. We used 
the RIDGE coherent network analysis algorithm for this study and show the sensitivity 
of this method to a variety of possible gravitational wave signals which could happen 
due to this pulsar glitch. Since the coherent method used here can calculate the 
detection statistics for the entire sky location, it is quite easy to search for glitch 
events from pulsars located at any sky position. We also studied the accuracy of the 
reconstruction of gravitational wave signals with RIDGE using matched filtering. The 
result of the matched filter analysis shows that the waveform reconstruction by RIDGE 
is nearly optimal, so our approach can provide information about the mechanism of 
pulsar glitches and a clue to infer the inner structure of neutron stars. 

Acknowledgments 

We would like to thank I Yakushin and Y Mochizuki for fruitful discussions and 
valuable comments on this paper. K.H. is supported by NASA grant NAG5-13396 to 
the Center for Gravitational Wave Astronomy at the University of Texas at Brownsville 
and NSF grant nsf-hrd0734800, SDM's work was supported by NSF grant PHY- 
0555842. S.Y. is supported by the Southeastern Louisiana University and NSF grant 
PHY 0653233. S.D. is supported by the Center for Gravitational Wave Physics at the 
Pennsylvania State University. The Center for Gravitational Wave Physics is funded 
by the National Science Foundation under cooperative agreement PHY 01-14375. T.S. 
is supported by a grant from the Office of Scholarly Research at Andrews University. 
This paper has been assigned LIGO Document Number P080032. 

References 



[1 

[2: 
[3; 

[4 

[6; 

[7; 

[s; 

[9; 

[10: 

[11 

[12: 

[is: 

[14 

[is: 
[16: 
[17: 
[is: 

[19 

[20: 



McCuUoch P, Hamilton P, McConnell D and King E 1990 Nature 346 

Baym G, Pethick C, Pines D and Ruderman M 1969 Nature 224 872 

Anderson P and Itoh N 1975 Nature 256 25 

Ruderman M, Zhu T and Chen K 1998 Astrophys. J. 492 267 

Takatsuka T and Tamagaki R 1988 Prog. Theor. Phys. 79 274 

Mochizuki Y, Izuyama T and Tanihata I 1999 Astrophys. J. 521 281 

Andorsson N and Comer G 2001 Phys. Rev. Lett. 87 

Itoh N 1983 Prog. Theor. Phys 69 3884 

Middleditch J, Marshall F, Wang Q, Gotthelf E and Zhang W 2006 Astrophys. J. 652 1531 

Abbott B et al 2004 Nucl lustrum. Meth. A 517 154 

Acernese F et al 2002 Class. Quantum Grav. 19 1421 

Flanagan C and Buchner S 2006 CBET 595 

Clark J, Heng S, Pitkin M and Woan G 2007 Phys. Rev. D 76 043003 

Tikhonov A N and Arsenin V Y 1977 Solutions of ill-posed problems (Washington, D.C.: V.H. 

Winston & Sons) 
Rakhmanov M 2006 Class. Quantum Grav. 23 S673 
Mohanty S 2002 Class. Quantum Grav. 19 1513 
Mukherjee S 2003 Class. Quantum Grav. 20 S925 

Hayama K, Mohanty S, Rakhmanov M and Desai S 2007 Class. Quantum Grav. 24 S681 
The VIRGO collaboration 2005 Class. Quantum Grav. 22 S1189 
Klimenko S, Mohanty S, Rakhmanov M and Mitselmakher G 2005 Phys. Rev. D 72 122002 



Searches for gravitational waves from pulsar glitches 9 

[21] Mohanty S D, Rakhinanov M, Klimenko S and Mitsebnakher G 2006 Class. Quantum Grav. 23 

4799 
[22] Barish B and Weiss R 1999 Phys. Today 52 44 
[23] Hayama K and Fujimoto M K 2006 Class. Quantum Grav. 23 S9 



matched filter output 



O 




I 



5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 



140 



matched filter output 



O 




I 







100 
fc 



120 



I 



7 



140 



